Shape fluctuations in the ground and excited + states 

of 30 Mg and 32 Mg 
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Large-amplitude collective dynamics ol shape phase transition in the low-lying states ol 30_36 Mg 
is investigated by solving the five-dimensional (5D) quadrupole collective Schrodinger equation. The 
collective masses and potentials of the 5D collective Hamiltonian are microscopically derived with use 
of the constrained Hartree-Fock-Bogoliubov plus local quasiparticle RPA method. Good agreement 
with the recent experimental data is obtained for the excited + states as well as the ground 
bands. For 30 Mg, the shape coexistence picture that the deformed excited + state coexists with 
the spherical ground state approximately holds. On the other hand, large-amplitude quadrupole- 
shape fluctuations dominate in both the ground and the excited + states in 32 Mg, so that the 
interpretation of 'coexisting spherical excited + state' based on the naive inversion picture of the 
spherical and deformed configurations does not hold. 
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Nuclei exhibit a variety of shapes in their ground and 
excited states. A remarkable feature of the quantum 
phase transition of a finite system is that the order pa- 
rameters (shape deformation parameters) always fluctu- 
ate and vary with the particle number. Especially, the 
large-amplitude shape fluctuations play a crucial role in 
transitional (critical) regions. Spectroscopic studies of 
low-lying excited states in transitional nuclei are of great 
interest to observe such unique features of the finite quan- 
tum systems. 

Low-lying states of neutron-rich nuclei around N = 
20 attract a great interest, as the spherical configura- 
tions associated with the magic number disappear in 
the ground states. In neutron-rich Mg isotopes, the in- 
crease of the excitation energy ratio E(A+)/E(2+) [J- 
H and the enhancement of B(E2; 2f -> of ) from 30 Mg 
to 34 Mg indicate a kind of quantum phase tran- 

sition from spherical to deformed shapes taking place 
around 32 Mg. These experiments stimulate microscopic 
investigations on quadrupole collective dynamics unique 
to this region of the nuclear chart with various theo- 
retical approaches; the shell model M - tlCI . the Hartree- 
Fock-Bogoliubov (HFB) method [ll|, [l2|], the parity- 
projected HF [IU, the quasiparticle RPA (QRPA) [lj, 
Il5| . the angular-momentum projected generator coordi- 
nate method (GCM) with [l6| and without [13, [3 re- 
striction to the axial symmetry, and the antisymmetrized 
molecular dynamics [191 ] - 

Quite recently, excited + states were found in 30 Mg 
[HSU and 32 Mg [H at 1.789 McV and 1.058 MeV, re- 
spectively. For 30 Mg, the excited Oj" state is interpreted 
as a prolately deformed state which coexists with the 
spherical ground state. For 32 Mg, from the observed pop- 
ulation of the excited 0j~ state in the (t,p) reaction on 
30 Mg, it is suggested [22| that the 0^" state is a spheri- 
cal state coexisting with the deformed ground state and 



that their relative energies are inverted at N = 20. How- 
ever, available shell-model and GCM calculations consid- 
erably overestimate its excitation energy (1.4 — 3.1 MeV) 
S E3, EE HH • It is therefore a challenge for modern mi- 
croscopic theories of nuclear structure to clarify the na- 
ture of the excited 0^" states. For understanding shape 
dynamics in low-lying collective excited states of Mg iso- 
topes near N = 20, it is certainly desirable to develop a 
theory capable of describing various situations in a uni- 
fied manner, including, at least, 1) an ideal shape co- 
existence limit where the wave function of an individual 
quantum state is well localized in the deformation space 
and 2) a transitional situation where the large- amplitude 
shape fluctuations dominate. 

In this article, we microscopically derive the five- 
dimensional (5D) quadrupole collective Hamiltonian us- 
ing the constrained Hartree-Fock-Bogoliubov (CHFB) 
plus local QRPA (LQRPA) method [23|. The 5D 
collective Hamiltonian takes into account all the five 
quadrupole degrees of freedom: the axial and triaxial 
quadrupole deformations (/?, 7) and the three Euler an- 
gles. This approach is suitable for our purpose of de- 
scribing a variety of quadrupole collective phenomena in 
a unified way. Another advantage is that the time-odd 
mean-field contributions are taken into account in eval- 
uating the vibrational and rotational inertial functions. 
In spite of their importance for correctly describing col- 
lective excited states, the time-odd contributions are ig- 
nored in the widely used Inglis-Bclyacv cranking formula 
for inertial functions. The CHFB + LQRPA method has 
been successfully applied to various large-amplitude col- 
lective dynamics including the oblate-prolate shape coex- 
istence phenomena in Se and Kr isotopes [23|, [24|, the 7- 
soft dynamics in sd-shell nuclei [25j , and the shape phase 
transition in neutron-rich Cr isotopes [26| . A preliminary 
version of this work was reported in Ref . (27j . 
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The 5D quadrupolc collective Hamiltonian is written 

as 

H C oii=T vib +T rot + F(/3,7), (1) 

T vib = \d pp (P, 7 )/3 2 + D^tf, 7 )/3 7 + ^ 77 (/3, 7 )7 2 , 

(2) 

1 3 

^--^^(AtK 2 , (3) 
k=l 

where T v ib and T rot are the vibrational and rotational 
kinetic energies, respectively, and V is the collective po- 
tential. The vibrational collective masses, Dpp, Dp y , and 
-D 77 , are the inertial functions for the 7) coordinates. 
The rotational moments of inertia associated with the 
three components of the rotational angular velocities 
are defined with respect to the principal axes. In the 
CHFB + LQRPA method, the collective potential is cal- 
culated with the CHFB equation with four constraints on 
the two quadrupole operators and the proton and neu- 
tron numbers. The inertial functions in the collective 
Hamiltonian are determined from the LQRPA normal 
modes locally defined for each CHFB state in the (/3, 7) 
plane. The equations to find the local normal modes 
are similar to the well-known QRPA equations, but the 
equations are solved on top of the non-equilibrium CHFB 
states. Two LQRPA solutions representing quadrupole 
shape motion are selected for the calculation of the vibra- 
tional inertial functions. After quantizing the collective 
Hamiltonian ([T]), we solve the 5D collective Schrodingcr 
equation and obtain collective wave functions 

*«i;ifG8,7,fi)= E $*iK(P,im\IMK), (4) 

K— even 

where <& a iK (/?, 7) are the vibrational wave functions and 
ifl\IMK) are the rotational wave functions defined in 
terms of V functions T>m K {Q). We then evaluate E2 
matrix elements. More details of this approach are given 
in Ref. 

We solve the CHFB + LQRPA equations employing, as 
a microscopic Hamiltonian, the pairing-plus-quadrupolc 
(P+Q) model including the quadrupole-pairing interac- 
tion. As an active model space, the two-major harmonic 
oscillator shells (sd and pf shells) are taken into account 
for both neutrons and protons. To determine the param- 
eters in the P+Q Hamiltonian, we first perform Skyrme- 
HFB calculations with the SkM* functional and the sur- 
face pairing functional using the HFBTHO code p8| . 
The pairing strength (Vq = —374 MeV fm -3 , with a cut- 
off quasiparticle energy of 60 MeV) is fixed so as to repro- 
duce the experimental neutron gap of 30 Ne (1.26 MeV). 
We then determine the parameters for each nucleus in 
the following way. The single particle energies are deter- 
mined by means of the constrained Skyrme-HFB calcu- 
lation at the spherical shape. The resulting single parti- 
cle energies (in the canonical basis) are then scaled with 
the effective mass of the SkM* functional m*/m = 0.79, 



since the P+Q model is designed to be used for single- 
particle states whose effective mass is equal to the bare 
nucleon mass. In 32 Mg, the N = 20 shell gap be- 
tween d 3 /2 and f?/2 is 3.7 MeV for the SkM* functional, 
and it becomes 2.9 MeV after the effective mass scal- 
ing. This value is appreciably smaller than the stan- 
dard modified oscillator value 4.5 MeV (29[. This spac- 
ing almost stays constant for 30 ~ 36 Mg. The strengths of 
the monopole-pairing interaction are determined to re- 
produce the pairing gaps obtained in the Skyrme-HFB 
calculations at the spherical shape. The strength of the 
quadrupole particle-hole interaction is determined to re- 
produce the magnitude of the axial quadrupole defor- 
mation p of the Skyrme-HFB minimum. The strengths 
of the quadrupole-pairing interaction are determined so 
as to fulfill the self-consistency condition [3(J. We use 
the quadrupole polarization charge Se po i — 0.5 for both 
neutrons and protons when evaluating E2 matrix ele- 
ments. We solve the CHFB + LQRPA equations at 
3600 /3-7 mesh points in the region < (i < /3 max and 
0° < 7 < 60°, with /3 max = 0.5 for 30 Mg and 0.6 for 

32,34,36 Mg _ 

Our theoretical framework is quite general and it 
can be used in conjunction with various Skyrme 
forces/modern density functionals going beyond the P+Q 
model. Then the effects of weakly bound neutrons and 
coupling to the continuum on the properties of the low- 
lying collective excitations, discussed in Refs. [13, EH, 
can be taken into account, for example, by solving the 
CHFB + LQRPA equations in the 3D coordinate mesh 
representation. However, it requires a large-scale calcu- 
lation with modern parallel processors and it remains as 
a challenging future subject. A step toward this goal has 
recently been carried out for axially symmetric cases [26| . 

Figure [T] shows the collective potentials V((3, 7) for 
3 Mg. It is clearly seen that prolate deformation 
grows with increase of the neutron number. The col- 
lective potential for 30 Mg is very soft with respect to (3. 
It has a minimum at (3 = 0.11 and a local minimum at 
/3 = 0.33. The barrier height between the two minima is 
only 0.24 MeV (measured from the lower minimum). In 
32 Mg, in addition to the prolate minimum at j3 = 0.33, 
a spherical local minimum (associated with the N = 20 
spherical shell gap) appears. The barrier height between 
the two minima is 1.0 MeV (measured from the lower 
minimum). The spherical local minimum disappears in 
34 Mg and 36 Mg, and the prolate minima become soft in 
the direction of triaxial deformation 7. In 34 Mg, the po- 
tential minimum is located at 7 = 10°. 

In Fig. [3J calculated excitation energies and E2 transi- 
tion strengths are compared with the experimental data. 
The lowering of the excitation energies of the 2f and 
states and the remarkable increase of B{E2\2\ — > 0^) 
from 30 Mg to 34 Mg are well described in this calcula- 
tion. The calculated ratio of the excitation energies, 
E(4f)/E(2f), increases as 2.37, 2.82, 3.26, and 3.26, 
while the ratio of the transition strengths, B(E2;4^ — > 
2+)/B(E2;2+ -> 0+), decreases as 2.03, 1.76, 1.43, and 
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FIG. 2: (Color online) Comparison of calculated excita- 
tion energies of the 2f and 4f states (upper panel) and 
B(E2;2+ -> Oi) values (lower panel) in 3 °- 36 Mg with ex- 
perimental data p]-Q. 



1.47, in going from 30 Mg to 36 Mg. Thus, the properties 
of the 2 + and 4 + states gradually change from vibrational 
to rotational with increasing neutron number. 

Let us next discuss the properties of the + states and 
the 2 + and 4 + states connected to the + states with 
strong E2 transitions. The result of calculation is pre- 
sented in Fig. [3] together with the recent experimental 
data. The calculated excitation energies of the + states 
are 1.353 and 0.986 MeV for 30 Mg and 32 Mg, respectively, 
in fair agreement with the experimental data [U, [22j . 
In particular, the very low excitation energy of the + 
state in 32 Mg is well reproduced. In our calculation, 
more than 90% (80%) of the collective wave functions 
for the yrast (excited) band members are composed of 
the K = component. Therefore we denote the ground 
band by 'the K = Oi band,' and the excited band by 'the 
K = O2 band.' The 2 + and 4 + states belonging to the 
K = O2 band appear as the second 2 + and 4 + states in 
30:32 Mg, while they appear as the third 2 + and 4 + states 
in 34 ' 36 Mg. Accordingly, we use 2 + 3 and 4 + 3 , to collcc- 



FIG. 3: (Color online) Excitation energies of the ex- 
cited O2 , and 4^3 states (upper panel) and the ra- 
tio B(£2;0 + -> 2+)/B(E2;0+ -> 2 + 3 ) of the inter-band 
E2 transition strengths between the K = O2 and K = Oi 
bands (lower panel). Experimental data are taken from 
Refs. [M[2l]. See texts for details. 



tively indicate the second or the third 2 + and 4 + states. 
The calculated ratios of the excitation energies relative to 
the excited 0^ state, [E(4 c + 3 )-E(0+ )}/[E(2+ 3 )-E(0+)], 
are 3.18, 2.87, 3.25, and 3.00, for 30 Mg, 32 Mg, 34 Mg, 
and 36 Mg, respectively. In the upper panel of Fig. [3] we 
also plot the rotor-model prediction for the excitation 
energies of the 4 + states estimated from the + — 2 + 
spacings in the K = O2 bands. The deviation from the 
rotor-model prediction is largest in 32 Mg indicating im- 
portance of shape-fluctuation effects. Although the cal- 
culated excitation spectrum of the K = 2 band in 30 Mg 
looks rotational, we find a significant deviation from the 
rotor-model prediction in the E2 transition properties. 
The calculated ratios of the E2 transition strengths, 
£(£2:4+3 -> 2+ 3 )/S(£2;2+ 3 -> 2 +), are 1.05, 1.54, 
1.47, and 1.51, for 30 ~ 36 Mg, respectively. The deviation 
from the rotor-model value (1.43) is largest in 30 Mg. The 
significant deviation from the simple rotor-model pat- 
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take into account not only the properties of the collective 
potential in the /3 direction but also its curvature in the 
7 direction and the collective kinetic energy (collective 
masses). This point will be discussed in our forthcoming 
full-length paper. As suggested from the behavior of the 
inter-band B(E2) ratio, the vibrational wave functions 
of the 2+ state are noticeably different from those of the 
0^ state in 30 Mg and 32 Mg, while they are similar in the 
case of 34 Mg. Next, let us examine the vibrational wave 
functions of the 0^ and 2^" 3 states in 30_34 Mg. It is im- 
mediately seen that they exhibit one node in the (3 direc- 
tion. This is their common feature. In 30 Mg and 32 Mg, 
one bump is seen in the spherical to weakly-deformed 
region, while the other bump is located in the prolately 
deformed region around f3 = 0.3 — 0.4. In 34 Mg, the node 
is located near the peak of the vibrational wave function 
of the 0^ state, suggesting that they have /3-vibrational 
properties. 




)3 0.6 f3 0.6 )3 0.6 



FIG. 4: (Color online) Vibrational wave functions squared 
2~2 K (A 7) 1 2 of the 0f,2+,0% and 2j 3 states in 

30_34 Mg. Contour lines are drawn at every eighth part of 
the maximum value. 



tern of the K = 2 bands in 30 Mg and 32 Mg, noticed 
above, can be seen more drastically in the inter-band 
E2 transition properties. In the lower panel of Fig. |3l 
we plot the ratio B(E2;0^ -> 2f)/B(E2;0f ->■ 2^ 3 ) of 
the inter-band transition strengths between the K = 0i 
and K = O2 bands. If the K = Oi and K = O2 bands 
are composed of only the K = component and the in- 
trinsic structures in the (/3, 7) plane are the same within 
the band members, this ratio should be one. These ra- 
tios for 34 Mg and 36 Mg arc close to one, indicating that 
the change of the intrinsic structure between the + and 
2 + states is small. In contrast, the ratios for 30 Mg and 
32 Mg are larger than 10, indicating a remarkable change 
in the shape-fluctuation properties between the + and 
2 + states belonging to the K = 0i and K = O2 bands. 

Figure 0] shows the vibrational wave functions squared 
J2k V&aiK (Pi 7)| 2 - Let us first examine the character 
change of the ground state from 30 Mg to 34 Mg. In 30 Mg, 
the vibrational wave function of the ground 0^ state is 
distributed around the spherical shape. In 32 Mg, it is re- 
markably extended to the prolately deformed region. In 
34 Mg, it is distributed around the prolate shape. From 
the behavior of the vibrational wave functions, one can 
conclude that shape fluctuation in the ground 0^ state is 
largest in 32 Mg. To understand the microscopic mecha- 
nism of this change from 30 Mg to 34 Mg, it is necessary to 
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FIG. 5: (Color online) (a) Vibrational wave functions 
squared, |$ ai j= ,if=o(/3, 7 = 0.5° )| 2 , of the 0j~ states in 
30_34 Mg. Their values along the 7 = 0.5° line are plotted 
as functions of /3. (b) Probability densities integrated over 7, 

P{0) ee rd 7 |<W =0 ,K=o(/3,7)| 2 |G(/3,7)| 1/2 , of the °t states 
in 30_34 Mg, plotted as functions of /3. (c) Same as (a) but for 

the O2 states, (d) Same as (b) but for the 0^" states. 

To further reveal the nature of the ground and excited 
+ states, it is important to examine not only their vibra- 
tional wave functions but also their probability density 
distributions. Since the 5D collective space is a curved 
space, the normalization condition for the vibrational 
wave functions is given by 

/ |1>a/K(/3,7)| 2 |G(/?,7)r /2 ^7 = 1 (5) 

with the volume element 

\G(fi,>y)\ 1 ' 2 d/3d>y = 2p 4 ^/W(/3, 7 )i?(/3, 7) sin 3 7 d/3d 7 , 

(6) 
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W{fi,i) ={DMP, l)D^{P,l) - [^ 7 (/3, 7 )] 2 }/3" 2 , 

(7) 

R(J3,j) =D X (/?, 7)^2 (fi, i)D z [fi, 7 ), (8) 

where -Dfc=i,2,3 are the rotational masses dehned through 
Jk = <i/3 2 D k sin 2 (7 - 2?rfc/3). Thus, the probability den- 
sity of taking a shape with specific values of (/3, 7) is given 
by Ek \® a iK(Pn)\ 2 \G{(i,l)\ 1/2 . Due to the /3 4 factor in 
the volume element, the spherical peak of the vibrational 
wave function disappears in the probability density dis- 
tribution. Accordingly, it will give us a picture quite 
different from that of the wave function. Needless to say, 
it is important to examine both aspects to understand 
the nature of individual quantum states. 

In Fig. [5J we display the probability density integrated 
over 7, P(/3) ee /d7|$ Q , /=0iK=0 (/3, 7 )| 2 |G(/3,7)| 1/2 , 
of finding a shape with a specific value of j3, to- 
gether with the vibrational wave functions squared 
\®a,i=o,K=o(fi, 7)| 2 for the ground and excited + states 
(a = 1 and 2). Let us first look at the upper panels 
for the ground states. We note that, as expected, the 
spherical peak of the vibrational wave function for 30 Mg 
in Fig. 03a) corresponds to the peak at f3 ~ 0.15 of the 
probability density in Fig. [SJb). In Fig. [5jb), the peak 
position moves toward a larger value of j3 in going from 
30 Mg to 34 Mg. The distribution for 32 Mg is much broader 
than those for 30 Mg and 34 Mg. 

Next, let us look at the lower panels in Fig. [5] for the 
excited states. In Fig. (5][c). the vibrational wave func- 
tions for 30 Mg and 32 Mg exhibit the maximum peak at 
the spherical shape. However, these peaks become small 
and are shifted to the region with /3 ~ 0.1 and f3 ~ 0.2 
in 30 Mg and 32 Mg, respectively, in Fig. [5jd). On the 
other hand, the second peaks at (3 ~ 0.3 and f3 ss 0.4 in 
30 Mg and 32 Mg, respectively, seen in Fig. [5{c) become 
the prominent peaks in Fig.[5jd). In 30 Mg, the bump at 
(3 ~ 0.1 is much smaller than the major bump around 
j3 ~ 0.3. In this sense, we can regard the Og state of 
30 Mg as a prolately deformed state. In the case of 32 Mg, 
the probability density exhibits a very broad distribution 



extending from the spherical to deformed regions up to 
(3 = 0.5 with a prominent peak at /? ~ 0.4 and a node 
at f3 ~ 0.3. The position of the node coincides with the 
peak of the probability density distribution of the the 
0^ state, as expected from the orthogonality condition. 
The range of the shape fluctuation of the 0^ state in 
j3 direction is almost the same as that of the 0^ state. 
Thus, the result of our calculation yields a physical pic- 
ture for the 0^ state in 32 Mg that is quite different from 
the 'spherical excited + state' interpretation based on 
the inversion picture of the spherical and deformed con- 
figurations. In 34 Mg, the peak is shifted to the region 
with a larger value of /3 and the tail toward the spherical 
shape almost disappears. 

In summary, we have investigated the large-amplitude 
collective dynamics in the low-lying states of 30 ~ 36 Mg by 
solving the 5D quadrupole collective Schrodinger equa- 
tion. The collective masses and potentials of the 5D 
collective Hamiltonian are microscopically derived with 
use of the CHFB + LQRPA method. Good agreement 
with the recent experimental data is obtained for the 
excited + states as well as the ground bands. For 
30 Mg, the shape coexistence picture that the deformed 
excited + state coexists with the spherical ground 
state approximately holds. On the other hand, large- 
amplitude quadrupole-shapc fluctuations dominate in 
both the ground and the excited + states in 32 Mg, so 
that the interpretation of 'deformed ground and spherical 
excited + states' based on the simple inversion picture of 
the spherical and deformed configurations does not hold. 
To test these theoretical predictions, experimental search 
for the distorted rotational bands built on the excited 0^ 
states in 30 Mg and 32 Mg is strongly desired. 
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